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ABSTRACT 

Nearby Gamma-Ray Bursts (GRBs) have been proposed as a possible cause of mass extinctions on 
Earth. Due to the higher event rate of GRBs at higher redshifts, it has been speculated that life as we 
know it may not survive above a certain redshift (e.g. z > 0.5). We examine the duty cycle of lethal 
(life-threatening) GRBs in the solar neighborhood, in the Sloan Digital Sky Survey (SDSS) galaxies 
and GRB host galaxies, with the dependence of the long GRB rate on star formation and metallicity 
properly taken into account. We find that the number of lethal GRBs attacking Earth within the past 
500 Myr epoch of the Ordovician mass extinction) is 0.93. The number of lethal GRBs hitting 
a certain planet increases with redshift, thanks to the increasing star formation rate and decreasing 
metallicity in high-z galaxies. Taking I per 500 Myr as a conservative duty cycle for life to survive, 
as evidenced by our existence, we find that there are still a good fraction of SDSS galaxies beyond 
z = 0.5 where the GRB rate at half-mass radius is lower than this value. We derive the fraction of 
such benign galaxies as a function of redshift through Monte Carlo simulations, and find that the 
fraction is ~ 50% at z ^ 1.5 and ^ 10% even at z ~ 3. The mass distribution of benign galaxies is 
dominated by Milky-Way-like ones, thanks to their commonness, relatively large mass, and low star 
formation rate. GRB host galaxies are among the most dangerous ones. 

Subject headings: Gamma Ray Burst 


1. INTRODUCTION 


Nearby high-energy transient sources have been con¬ 
sidered as a possible cause of mass extinction events on 
Earth as well as a potential threat to life in the fu- 


ture (IRudermanI 119741: iThorsettI 11995 

: IDar et alJ 119981: 

Gehrels et al.1120031: IMelott & Thomas 

120111: ISvensmarkI 

2012f). Among them, gamma-rav bursts IGRBsL the 


most violent explosions in the universe, have been re¬ 
garded as one of the most lethal high-energy transie nts 
to life (|Scalo &: Wheeled 120021 : iThomas et al.l[^05al lbll. 
With a 2-D simulation on how Earth’s atmosphere re- 
sponds to an intense 7-ray flux, IThomas et al.l (|2005ar) 
investigated the effect of a nearby (2 kpc away) GRB 
to life on Earth. A typical GRB at this distance has a 
7-ray fluence of 10® erg cm“^ (= 100 kJ m“^), which 
w ould cause seve r e dama ge to life on Earth. According 
to IThomas et al.1 (|2005all . such a GRB would lead to a 
rapid increase of nitrogen compounds (e.g. NO and NO2) 
in the atmosphere causing an on-average 35% ozone de¬ 
pletion in the stratosphere, increasing for years the solar 
UVB radiation flux at Earth. The resulting DNA dam¬ 
age, up to 16 times the normal level, is lethal to many 
organisms such as plankton, the base of the food chain. 
It may lead to extinction of creatures in higher trophic 
levels of the food chain due to starvation. Eurthermore, 
the opacity of NO2 in Earth’s atmosphere would result in 
a decreased surface temperature, which is speculated to 
be a ca use of a long-lasting ice age (jThomas et al.l[^n5al 
but see lThomas et al1l2015ll . The mass extinction in the 
late Ordovician (^ 447 M yr ago) could be due to the 
impact of a nearby GRB dMelott et al.l 120041 : iBambachl 
[2?M IMelott fc Thoma..sll2008ll . 

Long-duration G RBs are more dangero us than the 
short-duration ones (jPiran fc .limene^ 1201 41 1. Since long 
GRBs (LGRBs hereafter) are related to deaths of massive 


stars and the refore track the star formation history of 
the u niverse (iWooslev &: Blooiiil 120061: iKumar fc Zhand 
I2015I1 . it is naturally expected that GRBs become more 
lethal at higher red shifts where the s t ar for mation rate 
is higher. Indeed, iPiran fc .limenezl (j2014f l speculated 
that life as we know it cannot survive at z > 0.5 due to 
the frequent bombardment of GRBs at any location in a 
galaxy. 

One caveat in drawing such a conclusion is related to 
the duty cycle of lethal GRBs preventing development of 
advanced life. The time scale to re-develop advanced life 
is not well studied. Even though the Ordovician mass 
extinction happened ~ 447 Myr ago, more recent mass 
extinction events (e.g. the Cretaceous-Paleogene extinc¬ 
tion ^ 65 Myr ago that killed the dinosaurs) did not 
prevent re-emergence of advanced life forms (humans). 
Therefore it may be possible for life to tolerate a shorter 
(say, 50 Myr) lethal GRB duty cycle. Regardless of the 
biological details of how a GRB may kill life and how ad¬ 
vanced life forms re-emerge, our existence suggests that 
advanced civilizations can develop if lethal GRBs have a 
duty cycle comparable to the one inferred for the solar 
neighborhood in the Milky Way Galaxy. 

In this paper, we quantify the lethal GRB duty cycle 
in the solar neighborhood, and apply it as a conserva¬ 
tive life survival condition to study the “habitability” of 
various observed galaxies (e.g. Sloan Digital Sky Sur¬ 
vey [SDSS] galaxies and the GRB host galaxies). The 
aim is to address whether life can survive GRBs in the 
high-redshift universe. 

2. METHODOLOGY 

Eollowing IThomas et al.l l)2005all . we adopt Fc = 
10® erg cm“^ as a critical fluence defining a lethal GRB. 
Since GRBs trace star formation and most star forma- 
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tion happens in late type galaxies, we consider our test 
galaxies as disk galaxies with an exponential stellar mass 
column density where is the scale length. 

To the Hrst-order approximation, one may ignore varia¬ 
tion of specific star formation rate (sSFR) and metallicity 
within the galaxy. A GRB at tq from the galactic center 
defines a pair of cones where life is damaged. These cones 
have a radius R{L, Fc) = LAt/{A ttFc) and a solid angle 
n, where L is the peak isotropic luminosity of the GRB, 
and At is the rest-frame duration of the GRB, which is 
typically ^ 10 s. For a GRB at tq, the fraction of mass 
in the galaxy where life is damaged can be expressed as 


P (L, ro) =p[R{L, Fc), ro] = 




M* 


rQ + R 


y/R'^-ix-TQ)'^ 


dx 


E*e dy 


’•0-^ -yR-d-(a:-royd 


( 1 ) 


typical time scale T = 500 Myr is used to match ap¬ 
proximately the epoch of the Ord ovician Mass Extinc¬ 
tion of ^ (^447 - 44 3) Myr ago (|Sutchffe et al.l 120001 : 
IBrenchlev et al. I [2 00 311 . 

The parameter po denotes the average GRB event rate 
density in the local universe (z ~ 0). We know that the 
LGRB rate depends on star formation rate and metallic¬ 
ity. It is then relevant to introduce two correction fac¬ 
tors for the specific galaxy values with respect to the 
local average values. The specific SFR correction factor 
is introduced as 


sSFR 
“ sSFRo ’ 


( 4 ) 


where the l ocal specific star form ation rate is sSFRq = 
0.1 Gyr“^ (|Weinmann et al]l2011[l . The metallicity cor¬ 
rection factor is defined as 


4^([Fe/H]) 

Po([Fe/H]o)’ 


( 5 ) 


where /{, = |^ is the beaming correction factor. The 
fractionp(L, tq) is also the probability for a random GRB 
in the galaxy (including both those beaming towards and 
those beaming away from the planet) to kill life on a 
planet at rg from the galactic center. 

The high-luminosity LGRB luminosity function 
(LF), (f){L/L*) is characterized by a broken power 
law with a break luminosity L* ( Liaji g et al.l l2007t 
IWanderman fc Piranll20inHSun et "ni 201, 5h . e.g. 



The number of lethal GRBs that would attack the planet 
at To within a time duration T is 

F-i'max 

N= / f^^po^{L)p{L,ro)dL-TV{M,)UFRfFe,{3) 

^min 

where the range of integration is 10'’’® — lO^^erg s“^, the 
obseryed range of LGRBs, and po is the GRB event rate 
density above a specific luminosity in the local universe. 
Here the factor makes the correction from the ob¬ 
served GRBs (defined by the observed luminosity func¬ 
tion) to the total GRBs (including those not beaming 
towards us). Notice that it cancels out with the fb fac¬ 
tor in Eq. dm, so that the result does not depend on the 
poorly constrained parameter fb- 
Various studies reach a generally consistent conclu¬ 
sion, but with somewhat dif ferent parameters . Here we 
adopt the latest analysis by ISun et all l|2015ll using the 
largest long GRB sample. Since there is an evolution 
of the GRB luminosity function, we use the luminosity 
function and event rate density derived from the nearby 
(z < 1) sample with the following parameters (|Sun et al.l 
ai = 1.57, az = 1.8, L* = 4 x 10^^ erg s-\ and 
p{> 10®° erg s“^) = 1.6 Gpc“®yr“^. The occupied cos¬ 
mological volume of a galaxy with mass M* is estimated 
as V(M*) = M*/p*(z), where p*{z) is the average stellar 
density obtained by in tegrating the stellar mass function 
(|Mortlock et al.lF2015ll from 10^ — 10^^Mq for different 
redshifts, which is fit as p^,(z) = iQ^'^ ‘i6-o.39z Gpc“® 
(see also iMuzzin et al.l (|2013ll and papers therein). A 


Here P([Fe/H]) is the fraction of stars with metallic- 
ity poorer than [Fe/H] = —0.43 (equiv alent to 0.4 Zq) 
(|Bertelli et al.l[1994 IVirgili et al.l 1201 ill . assuming that 
the metallicity of the galaxy is a Gaussian distribution 
with the detected metallicity as the medium value and 
the st andard deviation aFe = 0.22, similar to the Milky 
Way (|Casagrande et al]l2011h . The function P([Fe/H]o) 
is a similar fraction, but with median value being the 
mean metallicity of the local universe [Fe/H]n = —0.006 
(|Gallazzi et al.ll2008l : IMadau &: Dickinsonll2014ll . 

The short GRB (hereafter SGRB) impact rate can be 
estimated with a similar method. The typical duration 
of a SGRB At is ^ 0.5 s. The luminosity function de¬ 
pends on the model of the merger delay time distribu- 
tion, but can be generally fitted as a single power law 
(|Sun et al.ll2015ll . We use the best-ht luminosity function 
parameters for a Gaussian delay time distribution, i.e. 
4>iL) oc with po(> 10®° erg s“^) = 1.3 Gpc“® yr“^. 
The range of integration is 10'’’° — 2 x 10®'’’ erg s“^ as 
observed. Since the SGRBs are likely due to mergers 
of two neutron stars or a neutron star-black hole sys¬ 
tem, the probability of a short GRB is no longer directly 
related to the fraction of the mass. Rather, we derive 
the probability by introducing an offset distribution of 
the afterglow with res pect to the host galaxies. Using 
the o bservational data (jFong et al.ll20ir)l: iFong fc Berged 
[Ml), we Hnd that the normalized (in units of the host 
galaxy scale length) offset distribution can be approxi¬ 
mated as a Gaussian function in logarithmic space, i.e., 
G(logic(r/rd), /i=0.32, (7=0.57). So the probability of a 
SGRB located in an area dA is 


p{L,ro) 


f G(logio(r/rd))dA 
(In 10)27rr2 


1 

/ G(logio(r/rd))d(logio(r/ra)) ’ 


( 6 ) 


where the integration in the numerator is over the dam¬ 
aging region around tq and the integration in the de¬ 
nominator is over the entire galaxy. Also there is no 
direct connection between the SGRB rate and sSFR or 
metallicity, so that /sSfr.sgrb = 1 and /pe.SGRB = 1 are 
adopted. 
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3. LETHAL GRB RATE IN MILKY WAY 



Fig. 1. — The number of lethal LGRBs (red solid line) and SGRBs 
(blue solid line) within 500 Myr as a function of distance vq from 
the Galactic Center. Two vertical black dashed lines indicate the 
half-mass radius (ri^/ 2 ) fhe Earth position, respectively. 


We first apply our methodology to the Milky Way 
(MW). The stellar mass of the MW is M:^ = 6.08 zb 
1.14 X 10^° M 0 (IMcMillanllMTI) . so that the occupied 
cosmological volume is V = 2.1 x 10~^ Gpc~^. With 
sSFR = 2.71 ± 0.59 x 10“^ Gyr“^ (iLicauia fc Newmaiil 
I2014fl . which is relatively small in the local universe, 
we get the sSFR correction factor /sSfr = 0.27. The 
metallicity distribution function of the MW is a Gaus- 
sian with the mean value [Fe/H] = —0.06 and crpe = 0.22 
l|Gasagrande et all 1201111 . it is relatively metal poor 
compared to the average value in the local universe, 
and the metallicity correction factor is /pe = 1-7. The 
half mass radius o f the MW is rd = 2.15 ± 0.14 kpc 
(jBovv fc Rixll201^ . and our Earth is locat ed at rn = 
8.33T0.35 kpc away from the Galactic center (|Reidlll99l : 
lEisenhauer h al.ll2003l iGillessen et al.l[200^ . 

The red line in Fig[T] denotes the number of lethal 
LGRBs for planets at different radii tq in the Milky 
Way within 500 Myr. At Earth’s location, the number of 
lethal LGRBs is 0.93 within 500 Myr. This is consistent 
with the hypothesis connecting the Ordovician mass ex¬ 
tinction with GRB activity. The number is larger than 
1 for To < 8 kpc. At the half mass radius, the number 
of lethal GRBs is Nm /2 = 2.91. The number increases 
to even larger values at smaller radii, suggesting that the 
regions close to the Galactic center are less habitable. 
Our existence suggests that N ^ 1 per 500 Myr can be 
regarded as a conservative criterion for the survival of 
advanced life forms. 

The blue line in Fig[T] shows the number of lethal 
SGRBs as a function of tq per 500 Myr. They are much 
rarer than LGRBs, suggesting that the LGRBs are the 
dominant species limit ing life in the universe (see also 
iPiran fc Jimenej 120141) . In the rest of the paper, we ig¬ 
nore the contributions of SGRBs. 


4. LETHAL GRB RATE IN OTHER GALAXIES 

Galaxies at high r edshifts have h igher sSFR 
(|Weinmann et al.1 120111: lllbert et M ] 120141 1 and lower 

metallicity (|Maiolino et al.l 120081 IMannucci et al.ll2010f) 
on average. These result in higher /pe and /sSfr 
correction factors, so that the number of lethal GRBs 
is expected to increase with redshift. In this section, 
we investigate the lethal GRB rate for other galaxies 
making use of the observational data directly. We 
employ the SPSS DR8 and DR1 2 galaxy samples 
(lAihara et al.ll201ll : iSmee et all 120131 : lAlam et aLll2015ll 
as the starting point, and extrapolate the sample to 
even higher redshifts with a Monte Garlo simulation, 
aiming at quantifying the redshift-dependence of the 
fraction of benign galaxies where life as we know it can 
survive GRBs. GRB host galaxies are also studied for 
comparison. For all the galaxies, instead of studying 
the ro-dependent lethal GRB rate, we only investigate 
the rate at the half-stellar-mass radius Vm /2 as a 
representation of the “habitability” of the galaxy. 


4.1. SDSS DR8 and SDSS/BOSS DR12 Samples 

The SDSS DR8 sample includes all spectroscopically 
classified galaxies in SDSS-l/B. It covers 9274 deg^ of 
the sky and includes 860,836 unique galaxies, with a me¬ 
dian redshift 0.12. The SDSS/BOSS DR12 sample is 
the SDSS-lB’s Baryon Oscillation Spectroscopic Survey 
(BOSS) galaxy catalog, which includes a LOWZ sample 
and a CMASS sample. There are 1,376,823 galaxies in 
total, with a median redshift 0.5. 

These samples are subject to observational selection 
effects. Since our primary aim is to check whether life 
can survive GRBs at high redshifts, for the first step we 
work on the observed sample only to make the case. 

We use the SDSS DR8 and SDSS/BOSS DR12 galaxy 
properties obtained from the Portsmouth GroupQ. In this 
sample, stellar mass, star formation rate, and metallicity 
are estimated by fitting the ugriz photo metric data with 
the p ublicly available code Hyper-Z (IBolzonella et al.l 
1200011. The stellar p o pulat ion models of IMarastonI ( 20051) 
and iMaraston et all (|2009ll are employed. Star formation 
rate is estimated with an exponentially declining star for¬ 
mation history (SFH), truncated SFH, or constant SFH 
model, whichever gi ves the best fit. The Initial Mass 
Function is Kroupa (|Kroupall200Ill . 

Figure [5] displays two 3-D maps of the galaxies in the 
samples. For clear illustration only galaxies with declina¬ 
tions in the range of [0°, 0.5°] are shown here. The left 
panel shows the galaxies with stellar mass 10^° Mq < 
M* < 10^^ Mq, and the right panel shows the galax¬ 
ies with 10^^ Mq < M* < 10^^ Mq. DR8 galaxies are 
shown as filled circles and DRI2 galaxies are shown as 
open circles. Golors are encoded with the estimated num¬ 
ber of lethal LGRBs at half mass radius N-m /2 in 500 
Myr. “Blue” and “green” galaxies have Ny [/2 < 1, so 
that they are safe for life. “Red” and “orange” galaxies 
have A^m /2 > 1, which are the most dangerous. “Yel¬ 
low” galaxies are MW-like galaxies. One can see that 
although “red” and “orange” galaxies gradually increase 
with redshift, there are still “blue” / “green” galaxies up 
to redshift 0.7. The naive speculation that no life can 


http: / / www.sdss.org/drl2/spectro/galaxy_portsmouth / 
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Fig. 2.— 3-D maps of SDSS DR8 (filled circles) and DR12 (open circles) galaxies with Dec in the range [0°, 0.5°]. Colors are encoded 
with the number of lethal GRBs within 500 Myr at the half mass radius each galaxy. The larger the AfM/ 2 i fhe more dangerous 

the galaxy. Left Panel: 10^*^ Mq < M* < 10^^ ^0; Right Panel: 10^^ Mq < M* < 10^^ ALq. 



Fig. 3. — The number of lethal GRBs within 500 Myr at the halt mass radius ia the galaxies at different redshifts. Blue points 

indicate the SDSS DR8 galaxies, and green dots indicate the SDSS/BOSS DR12 galaxies. Red stars are LGRB host galaxies and magenta 
squares are SGRB host galaxies. Left panel: 10^® Mq < M* < lO^t Mq; Right panel: lO^t Mg < M* < 10^^ ALq. GRB hosts are plotted 
in both panels regardless of their masses. 


survive beyond z = 0.5 is not confirmed. 

In order to take a closer look at how galaxies become 
more dangerous at higher redshifts, in Fig|3] we display 
Nm /2 of galaxies as a function of redshift. The left and 
right panels show the galaxies with stellar mass in the 
range of 10^° Mq < M* < 10^^ Mq, and 10^^ Mq < 
M* < 10^^ Mq, respectively. In each panel, DR8 galax¬ 
ies are shown as blue points and DR12 galaxies are shown 
as green points. The vertical dashed line shows the cor¬ 
responding redshift (z = 0.43) when Earth was formed. 
The horizontal solid line marks /Vm /2 = 1; above which 
most regions in the galaxy are dangerous for life. One 
can see in general for both mass ranges, galaxies become 
more dangerous at higher redshifts. There seems to be 
an offset (especially in the left panel) between the DR8 
(blue) and DR12 (green) galaxies, but it is caused by the 
different selection effects for the two samples. 


The fraction of galaxies with A^m /2 < 1 is shown in 
the left panel of Fig. SI The blue and green lines are 
for DR8 and DR12, respectively. In order to avoid the 
selection effects on luminosity and mass, we only focus 
on one mass bin < M^, < here. The 

results have similar trends for other mass bins. Each 
redshift bin is required to have more than 5,000 galaxies. 
Some galaxies are marked as SFR = 0. Most of them 
have a truncated SFH. Some others have a large age for 
an exponential decay of SFH as a function of time. For 
each redshift bin, the fraction of galaxies with SFR = 0 
are plotted as the dotted histograms. 

Since the DR8 sample tends to have more blue (late 
type) galaxies for a particular stellar mass bin whereas 
the DR12 sample originally searched for red (early type) 
galaxies, the two samples are both subject to a sample 
selection effect. In any case, the fraction lines (blue and 
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z z 

Fig. 4.— Left panel: The fraction of benign galaxies with A'^m/z < 1 as a function of redshift for the galaxies. SDSS DR8 galaxies 
(blue line) and SDSS/BOSS DR12 galaxies with 10^® < M < 10^^ Mq (green line) are shown. Monte Carlo simulated galaxies in the 
same mass bin (grey) are also shown for comparison. Right panel; The fraction of benign galaxies as a function of redshift according to 
Monte Carlo simulations. Different black lines denote different mass bins: from bottom up: 10® < M* < 10® Mq, 10® < M« < 10^® Mq, 
10^® < M* < 10^^ Mq and 10^^ < M* < 10^^ Mq, respectively. Red line shows the mass weighted fraction, which is similar to the result 
of mass bin 10^® < M* < 10^^ Mq. The shaded regions denote the uncertainty range of the sSFR factor between red and blue galaxies for 
each redshift bin (grey) or for mass-weighted case (orange). The green line shows a Gaussian fit with G(peak = 0.96, /r = —0.24, a = 1.44). 


green) derived from these two samples set the lower and 
upper limits of the fraction of benign galaxies as a func¬ 
tion of redshift. 

4.2. Monte Carlo Simulations 

In order to reduce the sample selection effects and in¬ 
vestigate redshift-dependence of the fraction of benign 
galaxies for life, we apply a Monte Carlo simulation 
aimed at achieving a more complete and unbiased sam¬ 
ple. 

We simulate blue and red galaxies separately since they 
follow different correla tions among redsh i ft, stellar mass, 
sSFR and metallicity (iKimm et al.ll200^ iMortlock et al.l 
1201511 . We simulate 100,000 blue galaxies in each Az = 
0.1 redshift bin. The number of red galaxies is es¬ 
timated by the ratio between the number of red and 
blue galaxies. The ratio is estimated by integrating 
the stellar mass function of blue and red galaxies from 
10^ — 10^^ M|q in each redshift bin (|Mortlock et al.l[2015l : 
iMuzzin et akl l^M^ . which gives an empirical relation 
logio(r/6) = -0.16-0.74^. 

For blue galaxies in each redshift bin, galaxy mass is 
given by stellar mass function as a function of redshift 
(jMortlock et al.ir2015[) . The sSFR distribution is consid¬ 
ered as a Gaussian distribution. The media n sSFR is a 
function of redshift and stellar mass (Eq.5 of lllbert et al.l 
l2014f) . The dispersion CTsSfr is given as 0.4 Gyr“^, a 
mean value of their Fig. 6. The metallicity distribution 
is also assumed to be a Gaussian function. The me- 
dian value is a functio n of stellar mass and SFR (Eq.4 of 
iMannucci et al.ll2010t ). In order to avoid the abnormal 
increase of very small and very large ^ 0 . 32 , the fourth- 
order term of that equation is dropped. A dispersion 
(jpe = 0.07 is used, which is the median value of the 
dispersions in those stellar mass-SER bins. 

The mass distribution of red galaxies as a function 


of redshift in IMortlock et al.l l)20I5ll is used to simulate 
the masses o f the red galaxies. T hey are mostly pas¬ 
sive galaxies (|Strateva et al.l[20Mll and their sSFRs are 
generally 0.01 to 0 .001 times of those of blue galaxies 
(jKimm et al.ll200'^ . One may assume that their sSFRs 
follow the same relation with stellar mass (jKimm et al.l 
HOO^, but have sSFRs a factor 0.001/0.003/0.010 of 
those of blue galaxies with the same masses. The ef¬ 
fect of different sSFRs of red galaxies are shown as the 
shadows in the right panel of FigH) 

The fraction of MG simulated galaxies with Nm /2 < 1 
is shown in the right panel of FigH] Black lines are for 
different stellar mass bins. From the bottom up, the lines 
are for 10® - 10^ Mq, 10^ - 10^° Mq, 10^° - lO^^ Mq 
and 10^^ — 10^^ Mq, respectively. The shaded regions 
enclose the results for different median sSFRs ratio be¬ 
tween red and blue galaxies. The center line is for 0.003, 
while the lower and upper ones are for 0.010 and 0.001, 
respectively. It reveals that the fraction of benign galax¬ 
ies in the small galaxy sample (M* < 10^° ATq) is low, 
which becomes smaller than 50% at z = 0.5. The frac¬ 
tion of benign galaxies increases with mass bins. For 
example, for M* > 10^° Mq, the fraction is still 50% 
at z; < 1.5. The red line and the orange shaded re¬ 
gion show the mass-weighted fraction of benign galaxies 
with Nyi /2 < 1- They are similar to the fraction line 
and shaded region for MW-like galaxies in the mass bin 
10^° Mq < M* < 10^^ Mq. This is because these galax¬ 
ies are massive and relatively common, and occupy most 
of the mass in each redshift bin. The fraction of benign 
galaxies is ~ 50% at z ~ 1.5, and is ~ 10% even at z ^ 3. 

4.3. GRB host galaxies 

LGRBs tend to reside in late type galaxies with 
high star formatio n rate and low metallicity (e.g. 
iSavaglio et al.l [200^ . whereas SGRB host galaxies tend 
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to be more diverse (|Fong &: Berged [20131 : lBergei][2014D . 
We study the GRB host galaxies, both for LGRBs and 
SGRBs, for their habitability. Table [T] lists the informa¬ 
tion of some LGRB and SGRB host galaxies with the 
desired information, including stellar mass, SFR, metal- 
licity, and scale length. Also listed are the fVM /2 values 
of those galaxies. They are also overplotted as red stars 
(LGRB host galaxies) and magenta squares (SGRB host 
galaxies) in Fig [3] One can see that LGRB host galaxies 
are among the extremely dangerous galaxies for life. The 
number of lethal LGRBs at half mass radius within 500 
Myr in all of the LGRB host galaxies is much greater 
than 1, suggesting that life can hardly survive in these 
galaxies. Meanwhile, SGRB hosts are relatively safer 
than LGRB hosts due to their much higher metallicity 
and lower sSFR. However, in all cases studied in our 
sample one has fVM /2 > 1, so that they still belong to 
dangerous galaxies. 

5. GONCLUSIONS AND DISCUSSION 

In this paper, we examined the duty cycle of lethal 
GRBs {Fc = 10® erg cm“^) in the Milky Way as well 
as observed galaxies at different redshifts. The duty cy¬ 
cle of lethal GRBs attacking Earth is about 1 per 500 
Myr, consistent with the time scale of the Ordovician 
mass extinction. Our existence suggests that such a duty 
cycle is long enough to allow advanced life (such as hu¬ 
man beings) to survive, so that such a duty cycle can 
be regarded as a conservative criterion for a benign envi¬ 
ronment. Adopting such a criterion, we investigated the 
fraction of benign galaxies as a function of redshift us¬ 
ing the SDSS DR8 and SDSS/BOSS DR12 samples and 
through MC simulations. We find that this fraction is as 
high as 99% in the local universe, suggesting that the cur¬ 
rent era is most suitable to the development of advanced 
life. As expected, the benign galaxy fraction decreases 
with increasing redshift due to the increase of SFR and 
decrease of metallicity at progressively higher redshifts. 
However, contrary to the naive expectation that no life 
can survive beyond z > 0.5, we show that benign galaxies 
do exist at redshifts much higher than 0.5. In particu¬ 
lar, the fraction is ~ 50% around z ^ 1.5, and is still 
^ 10% even at z ^ 3. In view that the birth of Earth 
itself corresponds to z 0.43, our results raise the excit¬ 
ing possibility that advanced civilizations are in principle 
not excluded even before the formation of the solar sys¬ 
tem, even though the fraction of “habitable” galaxies is 
smaller in the higher redshift universe. 

Finally, we discuss some possible uncertainties inher¬ 
ent to this analysis. 1. Many uncertainties are involved 
in making the argument that a GRB is lethal. These 
include how the atmosphere reacts to intense y-ray flux, 
how DNA/life reacts to exposure to UV light, and how 
alien life forms may be different from the ones of which 
we are aware. On the other hand, since our analysis fo¬ 
cuses on how life survives a GRB, these uncertainties do 
not affect our argument. The only biological connection 
used in this paper is our existence, suggesting that the 
local lethal LGRB event rate (1 per 500 Myr) can be 
regarded as a safe (conservative) duty cycle for life to 
survive. If the GRB damage is less severe, as sug gested 
by some recent studies (e.g. iThomas et al.l[2(HH ). or if 


advanced life can tolerate a higher lethal GRB rate (ar¬ 
gument presented in the second to last paragraph in in¬ 
troduction), then the fraction of benigh galaxies at high-z 
is even larger. 2. There exist uncertainties in the local 
GRB event rate density, which affect the estimated lethal 
GRB duty cycle at Earth. We used the latest results of 
iSun et all (|2015ll with the largest LGRB sample. Adopt¬ 
ing the GRB event rate results derived by other authors 
would give a local lethal GRB rate generally consistent 
with our value, with a difference at most a factor of 2- 
3. 3. The two galaxy samples included in our study 
are subject to substantial selection effects. The derived 
fractions, even with MC simulations, may not fully repre¬ 
sent the true fractions due to some un-modeled biases. In 
any case, we have argued that the DR8 and DRI2 sam¬ 
ples would bracket the true population. The fact that 
at least some galaxies are benign against lethal GRBs 
at high redshifts is robust, which is the main conclusion 
of this paper. 4. The dependence of LGRBs on SFR 
and metallicity is evidenced by the observations, but the 
exact correction factors (especially the metallicity one 
/pe) are not fully constrained. In any case, the general 
conclusion of this paper does not depend on the con¬ 
crete form of /sSFR and /pe. 5. We did not study the 
impact of SGRBs in other galaxies. For most star forma¬ 
tion galaxies, their contribution to lethal GRB number is 
negligible. However, for early-type galaxies where SFR 
is extremely low, SGRBs may dominate the lethal GRB 
attacking duty cycle. In any case, A^m /2 is always much 
less than unity in these early-type galaxies. For the pur¬ 
pose of this paper (to claim the fraction of benign galax¬ 
ies at high redshifts), ignoring the SGRB contribution 
does not affect the conclusions and is therefore justified. 
6. We did not study the life-damaging effect of other 
transients such as supernovae, which may have compara- 
ble damaging effects as GRBs in the loca l universe (e.g. 
iMelott fc Thom^l2011l : lSvensmarkll2012f l. While a more 
detailed study is needed to address the supernova effect, 
the general conclusion of this paper may not be modified 
when supernovae are considered. This is because super¬ 
novae also follow star formation history, so that they may 
be subject to the same /sSfr correction factor. On the 
other hand, their metallicity dependence is less apparent 
than LGRBs, so that their /pe correction factor, if any, 
should be smaller than that of LGRBs. Overall, their 
rate of increase with redshift is slower than LGRBs, so 
that their role in damaging life at high- 2 ; would be less 
significant than LGRBs. 
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TABLE 1 

GRB HOST GALAXY PROPERTIES 


GRB 

Z 

logM 

Mq 

sSFR 

Gyr-i 

12+ 

log(0/H) 

kpc 

'^M/2 

Ref 

970228 

0.695 

8.65 

1.19 

8.47 

2.458 

15.74 

1,2,3 

980425 

0.009 

8.68 

0.54 

8.16 

3.266 

5.14 

1,2,3 

980703 

0.966 

10.00 

1.66 

8.14 

1.342 

2696.83 

1,2,3 

990712 

0.434 

9.29 

1.23 

8.10 

1.591 

194.86 

1,2,3 

991208 

0.706 

8.53 

13.34 

8.02 

0.344 

3558.73 

1,2,3 

010921 

0.451 

9.69 

0.51 

8.15 

1.903 

157.67 

2,3 

11121 

0.362 

9.81 

0.35 

8.60 

12.671 

2.46 

2,3 

020903 

0.251 

8.87 

3.57 

8.22 

0.431 

927.26 

2,3 

030329 

0.168 

7.74 

2.00 

7.97 

0.804 

18.18 

2,3 

040924 

0.859 

9.20 

1.19 

8.23 

2.382 

117.08 

2,3 

SORB 

050709 

0.161 

8.80 

0.24 

8.50 

2.080 

3.08 

4,5 

051221A 

0.546 

9.40 

0.38 

8.80 

2.290 

3.62 

4,5 

061006 

0.438 

9.00 

0.24 

8.60 

3.220 

2.14 

4,5 

070724A 

0.457 

10.10 

0.20 

8.90 

3.640 

1.73 

4,5 

071227 

0.381 

10.40 

0.02 

8.50 

4.720 

4.75 

4,5 

130603B 

0.356 

9.70 

0.34 

8.70 

2.020 

14.36 

4,5,6 


References. — fl lBloom et al.| | |2002D ; ('2 lWainwright et al.| l |2007t) : 
(3 Savaglio et al.| ||2009D: f4 lLei^er & Berg^ | |2010D : fS lFong & Bergeij j2013f) : 
(6' et all l2^^) 
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